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ABSTRACT 


A detailed analysis of modelling and behaviour of buck boost converter has been done. The 
state space linear differential models are used to study the open loop time domain behaviour 
of the converter. To enable design of close loop feed back controllers continuous time 
representation of the converter is achieved by small signal linearized model which are 
obtained by linearizing the non linear state space averaged equations. Attempt is made to 
high light the fact that if switching information is neglected then state space average model 
can represent the converter behaviour for all kinds of perturbation. Using the small signal 
linearized model a state feed back controller is designed by pole placement and its response 
for perturbation in V re f, Vdc, R has been simulated. 
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CHAPTER I 


INTRODUCTION : 

In the field of industrial applications DC to DC converters play a vital role. There are many 
processes where direct human intervention is not possible or is not sufficient to meet the 
response time specifications that the plant under consideration demands. To implement a 
stable control loop for a system, it is required that system characteristics be known for both 
steady-state and transient. Therefore, detailed mathematical modelling of the system needs to 
be done using differential equations. 

The choice of the feedback variable to a great extent decides the close-loop response. 
However there is a method by which all the system variables that decide its state at any given 
time can be used to construct a control signal which constrains the system response to 
specifications, however the system must qualify necessary and sufficient conditions. The 
system under consideration may require more than one set of differential equations to model 
it completely, such that each set defines it during specific intervals of time. A continuous 
linear model for such a system has to be formed, so that either conventional or modem 
control techniques can then be applied to determine the close loop system parameters.. 



ORGANIZATION OF THE THESIS : 


The objective of the thesis was to model and study the characteristic of a Buck-boost DC to 
DC voltage converter and to implement a close-loop voltage control. Chapter 2 presents the 
study of buck-boost circuit, its operation under different modes of conduction and its 
mathematical model, both in differential equation and state space form. The state space 
model is then used to simulate open loop system responses for perturbation in duty ratio D 
and load resistance R. Chapter 3 explains the concepts and requirement of state-space 
averaging and how it is applied to the discrete differential equation models formulated in 
Chapter II, that define the converter in various modes of conduction to achieve continuous 
time representation of the converter. It is examined under what conditions do these models 
represent the actual converter or fail to do so. Also the linearization of continuous averaged 
models which are non-linear around the operating points so that design methods like Bode 
plot, Nyquist criteria or pole placement can then be used to determine the close loop 
parameters. In Chapter-IV how the pole location of the system governs the system dynamics 
and what conditions must the system satisfy such that its poles can be located at the desired 
location to <r. od/fv its dynamics, has been discussed. 


2 



CHAPTER n 


BUCK BOOST CONVERTER 

DC-DC voltage converters are used to obtain variable DC voltage from a fixed dc voltage 
source. In this chapter section 2.1 explains the buck-boost converter topology, circuit 
operation and how it can be represented by differential equations, Section 2.2 gives the 
state space representation of the converter in different modes of circuit operation and 
section 2.3 explains how the state space equations can be used to simulate converter 
operation in open loop. 

2.1 CONVERTER TOPOLOGY : - 

The work presented here is based on buck boost converter as shown in fig. 2.1. The 
switch S is turned on and off at frequency f s . In the discussion that follows the circuit 
operation is explained, equivalent circuits are formulated and differential equation model 
of the converter for different modes of operation is defined. 

MODE 1 : The switch S is on and the diode D e is reversed biased . The inductor L 
stores energy from the supply voltage V dc » the capacitor C discharges into load resistance 
R. The equivalent circuit is as shown in Fig. 2.2a. 

MODE 2 : The switch S is off and diode D e is forward biased. The equivalent circuit 
is as shown in Fig. 2.2b. The energy stored in inductor L in mode 1 is now transferred to 
capacitor C and load resistance R. In continous mode of conduction, mode 1 of the next 




MODE 2 : The switch S is off and diode D e is forward biased. The equivalent circuit is 
as shown in Fig. 2.2b. The energy stored in inductor L in mode 1 is now transferred to 
capacitor C and load resistance R. In continous mode of conduction, mode 1 of the next 
switching cycle starts before the inductor current in Mode 2 decays to zero. In discontinous 
mode of conduction there is a third mode. 

MODE 3 : In discontinuous mode of conduction the inductor current in mode 2 decays to 
zero before the next switching cycle starts. The inductor current tends to reverse but is 
blocked by diode D e . The equivalent circuit is as shown in figure 2.2c. Now both switch S 
and diode D e are off, inductor current remains zero till Mode 1 of next switching cycle starts 
and capacitor C continues to discharge in load resistance R. 

A special case of 2-mode Conduction is when the converter operates in 2-mode , but at the 
end of mode 2 the inductor current just becomes zero. This is called Boundary conduction. 
Figures 2.3 a, b, and c show the approximate waveforms of inductor current II and capacitor 
voltage V c during a switching cycle, along with initial and final values of II and V c for every 
mode during continous (2-mode) , boundary and discontinous (3-mode) operation of the 
converter respectively. 

As approximate nature of inductor current and capacitor voltage is known along with the 
equivalent circuits, converter operation can be represented by differential equations as 


follows : 





FIG 2.3(b): Approximate II &Vc Waveforms in Boundary Conduction Mode 



( 2 . 1 ) 


MODE 1 : Duration of this mode is to to ti or (DT C ) 



R,C 


Where Rj = R + r c 


1 = 



( 2 . 2 ) 


Solution to above equations are 


V c = V co e“ t/R|C 


V, V 

k = Tr+(ko-— )e- rtt,/L 


V L 


(2.3) 


(2.4) 


Where 0<t’<DT c . 

MODE 2: Duration of this mode is ti to t 2 or t f . Depending upon whether the cycle 

continous/discontinous It is assumed that load resistance is such that L-C circuit 
underdamped. 



-V e ri l 

Vc = - + k 


R,C R,C 


-V c r c /R, r L I L Rr c I L 

L L R, 


Solution to above equations is 


(C,Coscot"+C 2 Sincot") 


I L =e rt (C 3 Coscot"+C 4 Smcot" ) 


Where 0 < t” < (1-D)T C or Dj T c 


Where Ci = Vci 


r =1 Eu-v L+JL 

2 cd R,C C1 l R,C 


( 2 . 10 ) 


C 3 = Ili 


( 2 . 11 ) 
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( 2 . 12 ) 


C 4 =- 

CD 


M1 +1l Li3^ 

L L \ L r,u 


a = - 


1 

R,C 


+ ik + J^s_ 

L R,L 

2 


(2.13) 


co = 


4(R + r L ) r 1 r L Rr c V 
R,LC >,R,C L R,U 


2 


(2.14) 


MODE 3 : Duration of this mode is from t 2 to tf 


Vc 


R,C 


(2.15) 


II = 0 


(2.16) 


Solution to these equations is 
V c = V o2 e- , '" /RlC 
Il = 0 


(2.17) 

(2.18) 


0 <t”’ < (l-D-Di)Tc 



Given the values of V co and Ii 0 at begining of the switching cycle it can be determined 
whether present switching cycle for given duty ratio D will be continous / discontinous. 
(refer Appendix A-l for converter Parameters). If t 2 is the time taken by inductor current to 
decay to zero in mode 2 then depending upon its value the present switching cycle can be 
classified as a case of continous / discontinous/boundary conduction. 

Substituting Il 2 =0 in equation 2.8 (i.e. now value of t” = t 2 ). Results in 

e** 2 [C 3 Cosot 2 + C 4 Sincot 2 ] = 0 

since e 0 * 2 cannot be zero 

[C 3 Cosoot 2 + C 4 Sincot 2 ] = 0 

. S incot, C, 

Therefore 2. = — L_ 

Coscot 2 - C 4 

C 3 

tan cot, = — — 

-c 4 
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Therefore t 2 = 


tan~ 1 [C 3 /-C 4 ] 
co 


C4 could be positive or negative, to ensure that sign of C4 does not effect the value of t 2 



(2.19) 


If t 2 > (1-D)T C then the present switching cycle is continous. 

If t 2 = (l-D)Tc then the present switching cycle is a case of Boundary conduction. 
If t 2 < (1-D)T C then the switching cycle is discontinous. 

1 r 

, . , . 1 * 

2.2 STATE SPACE REPRESENTATION: wM 1 


A scalar differential equation can be transformed into state equations of the form 

X = Ax + Bu (refer Appendix A-2). Therefore equations 2.1, 2.2, 2.5, 2.6 and 2.15 and 2.16 
can be represented in the form 


X = A,X + B,u (2.20) 


X = A 2 X + B 2 u 


( 2 . 21 ) 



X = A3X + B3U 


( 2 . 22 ) 


Where X = 


V 


v 

c 

A. 

; x = 

c 

A. 


; u = v dc 


-1 

R,C 


; A2 - 


-1 

R,C 


R,C 



R 





R,C 


■ -1 

0 

-1 

r c R 


; A3 - 

R,C 


r, +- s — 



0 

0 

L 

L R.J 






B,= 


0 

LlJ 


;b 2 = 


;b 3 = 


Matrices Ai , A2, A3 , Bi, B2, B3 are called systems matrices. 


2.3 OPEN LOOP TIME DOMAIN SIMULATION OF THE CONVERTER : 


Using the equations 2.20 to 2.22 converter operation in continous or discontinous 
mode can be simulated. As seen the value of t2 decides whether present switching cycle is 
continous/discontinous. Algorithm used for simulation is shown in figure 2 . 4 . In open-loop 
converter operating point is defined by Duty Ratio D, input voltage Vd C and load resistance R. 



Define Initial Operating Point 
&Record Initial Conditions(Vco 
&Ilo) 



FIG 2.4 : Open - loop Time Domain Simulation Using State - Space Differential Ean Model 










Cases given in Table 1 have been considered for converter simulation in open loop. 


S.No. 

Initial Condition 

Final condition 

■1 

Duty Ratio 

Load Resistance 

Duty Ratio 

Load 

Resistance 

1 

0.1 

10ft 

0.1 

5ft 

2 

0.1 

10ft 

0.1 

20ft 

3 

0.5 

10ft 

0.5 

30ft 


0.1 

30ft 

0.3 

30ft 

5 

0.1 

5ft 

0.01 

5ft 

6 

0.1 

10ft 

0.1 

30ft 

7 

0.1 

30ft 

0.1 

50ft 

8 

0.1 

30ft 

0.15 

30ft 


Table 1 

Results of simulation are shown in Fig. 2.5 to Fig. 2.12. Every operating point defined by D, 
V dc and R has a unique value of capacitor voltage and inductor current in steady state. If the 
ripple is neglected the average capacitor /output voltage and inductor current can determined 
without solving the differential equation (discussed in chapter 3). From the time domain 
simulation of the converter in open-loop maximum overshoot , delay time, rise time and 
settling time for perturbation in D and R can be observed. Also the ripple in inductor current 
and output voltage can be observed during transient and steady state. Result show that as the 
load resistance is increased the converter tends to discontinous operation during transient and 
steady state, and depending upon its value as the value of duty ratio D is decreased the 
converter either goes into discontinous operation or remains in continous operation, and as 
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FIG 2.5: System Response for Change in R from 10 to 5 Ohms for D=0.1 
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the value of D is increased the converter either remains in discontinous operation or tends to 
continous operation. Hence it can be concluded that zones of continous, discontinous and 
Boundary operation can be defined, (discussed in chapter 3). 



CHAPTER m 


STATE SPACE AVERAGE MODEL 

The characteristics of switched systems can be studied by averaged models using simpler 
programming techniques that give faster output compared to the linear differential equation 
model. However certain information is obscured by such models. Section 3.1 discusses the 
state space averaged model representation of the buck-boost converter in the 2-mode and 3- 
mode operation. Section 3.2 uses the state-space averaged model to highlight characteristics 
of the buck-boost converter. Section 3.3 discusses how state-space averaged models 
represent the converter operation for all kinds of perturbations in Duty ratio D and load 
resistance R. Since state-space averaged models are non-linear, they have to be linearized 
about a DC operating point, so that feed back design techniques can be used to provide close- 
loop control for the converter. Section 3.4 discusses how small signal approximation gives 
linearized model of the converter about a DC operating point. 

3.1 STATE SPACE AVERAGED MODEL : 


To design a closed-loop system for the converter for voltage regulation the linear 
differential equation model is of little help. Since none of the equations 2.20, 2.21 or 2.22 
describe the converter for time interval long enough to include a discontinuity. State space 



averaging allows discontinuous (switched) systems to be approximated as a continuous but 
non-linear system (refer Appendix B-l). Linearization allows the resulting non-linear system 
to be approximated as a continuous and a linear system and the resulting equations can then 
be used to determine parameters of close-loop feedback system for desired transient and 
steady state response (refer appendix B-2). The continuous time state space averaged 
representation of Buck-boost converter is given by 


X = AOX + BOu (3.1) 

Where A0=Ai D + A 2 (1-D) (3.2) 

-1 R(1 - D) 

R,C R,C 

-R(l-D) -r L -r c R(l-D) 

R,C L LR, 

Where B0=Bi D + B 2 (1-D) (3.3) 

"O’ 

= D 

.L. 


During two mode or continuous mode of conduction, and 


A0 = AiD + A 2 D 1 +A 3 (1-D-D 1 ) (3.4) 

J± R^ 

R,C R,C 

-RDi a(D + Di) , r c RD, 

R,L ^ L LR, 


BO = Bi D + B 2 Di + B 3 (1-D-Di) 


(3.5) 



„ n _ Duration of mod e2 

Where Di = 

T 

C 


during three-mode or discontinuous mode of conduction. 


In feed back control systems D is a function of X and u. Thus continuos approximation of 
two or three switched linear systems is a non-linear system. Another advantage of using 
continuous time state space averaged model is that behaviors of the converter in time-domain 
can be studied (neglecting the switching information that can be studied exclusively by using 
only differential equation model) by using simpler programming techniques that give faster 
output compared to complete time domain analysis for perturbation in duty ratio D and load 


resistance R. Now X= 



represents average value of state variables. 


3.2 CHARACTERISTICS OF BUCK-BOOST CONVERTER : 


As observed in chapter-II the converter operation can be distinctly classified into 
zones of continuous, boundary and discontinuous operation, depending upon the value of 
load resistance R and duty ratio D. Also it was observed that if load resistance is kept 
constant and value of duty ratio D is increased the converter tends to boundary and then 



continuous mode of conduction, or if duty ratio D is kept constant and load resistance R is 
increased then the converter transits from continuous mode of conduction through boundary 
to discontinuous operation. In this section for given load resistance R how the classification 
of the three zones can be done on basis of duty ratio D, also in the discontinuous zone, for 
given value of R and D, how is Dj effected and how for given load resistance R the average 
capacitor voltage and inductor current vary with duty ratio D, is examined. 

3.2.1 VALUE OF DUTY RATIO D AT BOUNDARY OPERATION (D houaistry ) : 

(Asa function of load resistance R) 

During the two mode operation of the converter for given load resistance R there is a 
unique value of duty ratio D for which the inductor current at end of mode-2 just becomes 
zero. In steady state. 


0 



therefore equation 3.1 becomes. 

0 = AOX + BOu 

Substituting the values of AO and BO from equations 3.2 and 3.3 since boundary conduction s 
a special case of two mode conduction, results in two simultaneous equations : 


-1 

RC 


V e + 


R(l-D) 

RjC 


k 


= 0 


9PS 



and 


-R(i-D) 

R,L 


V. 


r L , r c R(l-D) 


— + - 
L 


R,L 


'^r v -o 


solving these equations gives 


V = 


DRR.Vdc 


R 2 (l-D) + r c R + 


(1-D) 


3(A) 


Where V c = Average Capacitor Voltage 
= V 0 (Average output Voltage) 

Consider a simple tl -L Ckt as shown in Figure 3.1 and consider that initial value of the 
inductor current at t=0 is iL =0 . This represents the equivalent circuit of the converter in 
Mode 1 of boundary operation at the input side. 

Then value of ii, at time t>0 is given by 

i L =— (l-e' rLt/L ) 
r L 

Since lim e“ x = 1 - x 

x-»0 

■l = 

r L -»-o r L L 

Using the above relation 
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In = Peak inductor current at end of Mode 1 


= V <, c A DT ° 

r L L 

= V do DT c 

L 


(3.6a) 


From figure 2.3 (b) it is observed that average value of inductor current during a switching 
cycle in Boundary conduction is 


Substituting value of In gives 


Ilb - 


1 V dc DT c 

2 L 


(3.6b) 


Equating the input-output power for Boundary conduction 
VdclLB DT C = V 0 Io T c 
Where I 0 = V 0 /R 

substituting the values of Ilb , Io and V 0 gives a fourth order quadratic equation. 

RV + 2R 3 r c u 3 + (2R 2 R,r L + R 2 r c 2 - 2Lf s RR 2 )u 2 + 2RR,r c r L u + R 2 r 2 = 0 (3.7) 
Where u = (l-D) 

One of the roots of above equation gives the value of ( 1-Dboundaiy) for given load resistance 
R, in steady state. 

Therefore D b0 undaiy = 1-u 

If in equation 3.7 r c = tl =0 then 



□boundary 


O' 


Vdc 



n 


L 


a 


FIG 3.1: Simple it - L Circuit 



FIG 3.2: Dboundary Versus Load Resistance R 




R 4 u 4 + (-2Lf s R 3 )u 2 = 0 


Therefore Dboundary 1 

which implies that Dboundary is proportional to L/R ratio. Also if Dboundary is made zero 
(implies only continuous operation of the converter is possible for given load resistance R) 
then the value of load resistance R above which both continuous / discontinuous operation of 
the converter is possible is directly proportional to value of L and £ . 

Consider the following cases : 

(a) A plot of R versus D boU ndary for L=100pH and f s =100KHZ 

(b) A plot of R versus Dboundary for L=10pH and f s = 100KHZ 

(c) A plot of R versus Dboundary for L=10pH and f s = 100MHZ 

(d) A plot of R versus Dboundary for L=1000pH and f s =100KHZ 

As shown in Figures 3.2 

It is observed that for some values of load resistance R, value of Dboundary is negative, it can be 
concluded that for these values of load resistance the converter shall never operate in 
discontinuous mode in steady state, however low is the value of D. Larger the values of 
inductance L the converter remains in continuous operation for given value of D for larger 
values of R in steady state. Also Dboundary remains same for same L/R ratio for different 
values of L. The value of R above which both continuous and discontinuous operation of 
converter is possible is proportional to L and f s . 




For given value of L and Vd C the area below the Dboundary curve represents all operating points 
that define the discontinuous operation zone of the converter in steady state and area above 
the curve defines all operating points lying in the continuous operation zone of the converter 
in steady state. 

3.2.2 DURATION OF MODE- 2 : ( D! T c ) 

Figure 2.3(c) shows that during discontinuous conduction the inductor current decays 
to zero and remains zero till the beginning of the next switching cycle, but in continuous 
conduction since the next switching cycle starts before the inductor current in Mode-2 decays 
to zero therefore duration of Mode-2 is taken as (1-D)T C and for boundary condition as (1- 
Dboundary) T c . However it is to be seen that for load resistance R for which both continuous 
and discontinuous operation is possible, how is duration of Mode-2 and Mode-3 effected as 
the value of D is increased. From figure 2.3(c) the capacitor voltage at end of Mode-2, V C 2 
can be expressed as : 

Vc2 = V cl e (1 ' Dl)T ‘ /T ' c (3.8) 


Substituting this value of Vc 2 in equation 2.7 and t” = DiT c results in 
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(3.9) 


V cl e (1-D ' )T/R,c = e CTD ' T ‘ [C 3 CoscoD,T c + C 4 Sin©D,T c ] 

In equation 2.8 if I L2 = 0 and t” =DiT c then 

e aD ‘ Tc {C 3 CoscoD,T c + C 4 SincoD,T c } = 0 

since e aD,Tc cannot be zero. 

C 3 CoscoD,T c + C 4 SincoD,T c = 0 (3.10) 

Substituting the values of C 3 and C 4 in equation 3.10 results in 


V c , (l-r./RJSincoD.T, 

C ° L | CoscoD.T, -(o + f + ||-)Sin©D 1 T c 
I JL KjL 


CO 

Substituting the values of C 3 , C 4 and Ili in equation 3.9 and eliminating V c i results in 


e (l-D,)T c /RiC _ e aT c D, 


Cos(coD,T c ) + 
R(l-r c /R.) 


SinV^D, 


R,Coo L 


, h. R^c \ 

(a + J, + ^L) 


Cos((»D 1 T c ) - 


L Ri r c 


CO 


-SmCfflDjTJ 


(cm- ) 

v p c 

— Sin(coD,T c ) 


co 


= 0 


(3.11) 
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FIG 3.3:(a) Duration of Mode2 as given by eqn(3.1 1) 




For given parameters of the converter and load resistance R solution of the 
transcedental equation 3.11 will give us the value of Di and DiT c (the duration of Mode-2) or 
the time required for inductor current in Mode-2 to decay to zero, but in continuous mode of 
conduction duration of Mode-2 is (1-D)T C . As seen in Figures 3.3 (a & b) for the following 
cases. 

(a) Plot of D VsDi forR = 20Q 

(b) Plot of D VsDi for R = 30 £2 

(c) Plot of D VsDi forR = 40Q 

(d) Plot of D VsDi forR = 50Q 

in Figure 3.3(a) to Figure 3.3(b) 

It is observed that value of Di remains constant for given value of R, irrespective of 
values of D. This is so because as seen from equivalent ckt in Figure 2.2(b) the rate of decay 
of inductor current is decided by the value of R and capacitor C. But in actual operation for 
D > Dboundary, Di= (1-D). The results show that in discontinuous zone Di=(l-Db 0U ndary) where 
Dboundary for given load resistance R is given by Equation 3.11. Therefore it can be concluded 
that in steady state during the discontinuous operation for given value of R it is only the 
duration of Mode-1 and Mode-3 that change while that of Mode-2 remains constant with 
change in value of D. 
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With this information given any operating point a state space averaged model of the 
converter can be defined, and the variation in steady state average values of Inductor current 
and Capacitor voltage for given load resistance and can be studied, as shown in Figures 
3.4(a) and (b). Simultaneous equations in section 3.2.1 give the steady state values of V c and 
II during continuous mode of conduction, similarly we can obtain relations for discontinuous 
mode of conduction by substituting in 0=A0X+B0u values of AO and BO as given by 
equations 3.4 and 3.5 as follows. 


— V C +^I L =0 
RjC R,C L 


-RAV d c 

R,L 


L (D + Di ) + ^ 


RjL 


I L 


-DVdc 

L 


Solution to the above simultaneous equations gives 


V = V = 


RR,DV dc 


R'D,+R,r L (l + ^-) + r t R 


(3B) 


Else using the relation 0 = AOX + BOu 

the steady state value of the average capacitor/output voltage and inductor current 
corresponding to given D and R can be determined as 


X0 = 



= inv(A0)B0V do 
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AVERAGE INDUCTOR CURRENT(IL) AVERAGE OUTPUT/CAPACITOR VOLTAGE(Vc) 



FIG 3.4:(a):Average Output Voltage versus Duty Ratio - '' />; 





Where AO and BO are given by equations 3.2 and 3.3 or 3.4 and 3.5 depending upon the 
operating point (defined by D,R and Vdc ) lies in continuous or discontinuous operating zone 
respectively. 

It is observed that for a given load resistance the output voltage and average inductor current 
increase gradually and marginally as D is increased, till D=0.8. Also for different values of 
load resistance R the output voltage Vo is not significantly different for given value of D till 
D=0.8. The output voltage varies significantly for different values of the load resistance R in 
the range D=0.8 to 0.99. For all values of R the slope of II vs D curve increases sharply 
beyond D=0.8. 

3.3 CONVERTER SIMULATION USING STATE-SPACE AVERAGED MODEL : 

If the converter behaviour has to studied in time domain for perturbation in duty ratio 
D and load resistance R without observing the switching information than state space average 
model provide means of doing so by employing simpler programming techniques which also 
provide faster output. Consider the cases give in Table- 1. As observed in Chapter-II 
converter operation can be classified as one of the following cases. 
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SI. No. 

Initial condition 

Transient State 

Final Condition 

1 . 

Continuous Mode 

Continuous Mode 

Continuous Mode 

2. 

Continuous Mode 

Discontinuous Mode 

Continuous Mode 

3. 

Continuous Mode 

Discontinuous Mode 

Discontinuous Mode 

4. 

Discontinuous Mode 

Discontinuous Mode 

Discontinuous Mode 

5. 

Discontinuous Mode 

Continuous Mode 

! 

Discontinuous Mode 

6. 

Discontinuous Mode 

Continuous Mode 

Continuous Mode 


Consider case 1 of Table- 1, if we use only the continuous mode state space averaged 
model to simulate the perturbation, we observe that if the ripple is neglected the results 
match those of the differential equation model very accurately, as shown in the Fig. 3.5 (a & 
b) . However, if the same model is used to simulate the case 5 we observe that the results are 
incorrect as shown in Fig. 3.5 (c & d). In the differential equation model the value of t 2 
decides whether the present switching cycle is continuous or discontinuous, but in state space 
averaged model t 2 holds no validity. Therefore, the conditions that are to be examined while 
using state space averaged model to effect the model transition, so that the model accurately 
represents the process for all cases of perturbations in load and duty ratio are the value of II 
and Dj. Equation 3.6(b) gives the value of average inductor current Ilb for a given value of 
duty ratio D when the converter is operating at the boundary condition. A value of average 
inductor current I L below this value at any instant for given D implies the converter may be in 
discontinuous conduction. The converter operation while using averaged model can be 
defined in Il-Di plane as shown in Fig. 3.6. 































FIG 3.5(a)! System Response for change in R from 10 to 
5 Ohms for D = 0.1;using state - space averaged model 
defined by eqns (3.2) & (3.3) 


FIG 3.5(b)! System Response for change in R from 10 
to5 Ohms for D = 0.1;using state - space averaged model 
defined by eqns (3.2) & (3.3) 







FIG 3.5(c)! System Response for change in D from 0.1 
to 0.01 for R = 5 Ohms ; using state - space averaged 
model defined by eqns (3.2) & (3.3) ; results show Model 
failure 


FIG 3.5(d)! System Response for change in D from 0.1 
to 0.01 for R = 5 Ohms ; using state - space averaged 
model defined by eqns (3.2) & (3.3) ; results show Model 
failure 


FIG 3.5 







if in equation (3 B), r c = r L = 0 then 


V c = 


VdcD 

D, 


Therefore 


D,= 


V dc D 

V c 


(3.12A) 


If the output of an averaged model at any instant satisfies the conditions : 


II < Ilb 
D,<(l-D) 

then for further solution averaged model as defined by equations 3.4 and 3.5 is used till the 
above conditions are satisfied. 

Case 2 to 5 given in Table 1 are simulated, using the above logic for Model transition. 
The results are shown in figures 3.7 to 3.10 respectively. 

It is observed from Fig. 3.10 (b) that despite implying model transition conditions the 
model fails because the above conditions of Model transition are not applied to the 
intermediate output of ode'45' subroutine. Fig. 3.10 (c & d) shows the response for case 5 of 
Table 1 when conditions of model transition are checked at every intermediate output of 
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FIG 3.6: II- Dl Plane: Conditions of Model transition fromCM-O- DCM 
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FIG 3.7(a)! System Response for change in R from 10 
to 20 Ohms for D = 0.1 ; using conditions given by 
eqns(3.6 b) and (3.12 A)for model transition. 



FIG 3.7(b) : System Response for change in R from 1 0 
to 20 Ohms for D = 0.1 ; using conditions given by 
eqns(3.6 b) and (3.12 A)for model transition. 
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FIG 3.8(a); System Response for change in R from 10 
to 30 Ohms for D = 0.5 ; using conditions given by 
eqns(3.6 b) and (3.12 A)for model transition. 


FIG 3.8(b)! System Response for change in R from 10 
to 30 Ohms for D = 0.5 ; using conditions given by 
eqns(3.6 b) and (3.12 A)for model transition. 
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FIG 3.9(a): System Response for change in D from 0.1 FIG 3.9(b): System Response for change in D from 0.1 

to 0.3 for R = 30 ; n si ng conditions given by eqns(3.6 b) to 0.3 for R = 30 ; using conditions given by eqns(3.6 b) 

and (3.12 A)for model transition. and (3.12 A)for model transition. 




ode45 subroutine. Figures 3.1 1 to 3.13 show the results for cases 6, 7 and 8 given in Table 
1, when the model transition conditions are implied to every intermediate output of the 
‘ode45’ subroutine. It is observed that the state space average model fails to represents the 
transient process, when the results are compared with those of differential equation model. 

Consider the relation given by equation (3B) 

RR,DV dc 

R 2 D| +R,r L (l +^-) + r c R 

Solving for Di from the above relation results in a quadratic equation. 


R 2 D, 


( + R,r L + Rr c Id, + R,r L D = 0 


(3.12B) 


One of the roots of the equation (3.12B) gives the value of Di for given R and D and 
Capacitor voltage V e . Figure 3.14 to 3.16 show the result for cases 6, 7 and 8 of Table 1 
when model transition condition are implied using the values of Ilb and Di as given by 
equations (3.6b) and (3.12B) to the intermediate output of ‘ode45’ subroutine. The results 
show that the state space average model is very sensitive to conditions implied for model 
transition. Figure 3.17 shows that for case 5 in Table 1 while the value of Di given by 
equation 3.6 (a) gives successful results, the value of Di given by equation 3.6 (b) represents 
the current transient incorrectly. Therefore apart from the value of Di what are the other 
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FIG 3.10 (a) 


FIG 3.10 (b) 


FIG 3.10:System Response for change in Duty Ratio D from 0.1 to 0.01 for 
R = 5 Ohms using Model transition criteion defined by eqns(3.6b & 3.12A) 




FIG 3.10(c) FIG 3.10(d) 


FIG 3.10:System Response for change in Duty Ratio D from 0.1 to 0.01 for 

R = 5 Ohms using Model transition conditions defined by eqns(3.6b & 3.12A) and 
implying it to intermediate output of ‘ode45’ subroutine 






CAPACITOR VOLTAGE(Vc) 



FIG 3.11(a) 



FIG 3.11(b) 


FIG 3.11 :System Response for Change in R from 10 to 30 for D = 0.1 ;nsing Model 
transition conditions as defined by eqns(3.6b) & (3.12A) and implying to 
intermediate output of ‘ode45’ subroutine. 




FIG 3.12(a) FIG 3.12(b) 

FIG 3.12 : System Response for Change in R from 30 to 50 Ohms for D = 0.1 ;using 
Model transition conditions as defined by eqns(3.6b) & (3. 12A) and implying to 
intermediate output of ‘ode45’ subroutine. 





CAPACITOR VOLTAGE(Vc) 




FIG 3.13(b) 


FIG 3.13 :System Response for Change in D from 0.1 to 0.15 for R = 30Ohms;using 
Model transition conditions as defined by eqns(3.6b) & (3.12A) and implying to 
intermediate output of ‘ode45’ subroutine. 




FIG 3.14 rSystem Response for Change in R from 10 to 30 for D = 0.1 ;using Model 
transition conditions as defined by eqns(3.6b) & (3.12B) and implying to 
intermediate output of ‘ode45’ subroutine. 
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FIG 3.15(a) FIG 3.15(b) 

FIG 3.15 : System Response for Change in R from 30 to 50 for D = 0.1 ;using Model 
transition conditions as defined by eqns(3.6b) & (3.12B) and implying to 
intermediate output of ‘ode45’ subroutine. 



FIG 3.16(a) FIG 3.16(b) 


FIG 3.16 -.System Response for Change in D from 0. 1 to 0. 15 for R = 30 Ohms using 
Model transition conditions as defined by eqns(3.6b) & (3.12B) and implying to 
intermediate output of ‘ode45’ subroutine. 
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FIG 3.17(a) FIG 3.17(b) 

FIG 3.17 : System Response for Change in D from 0.1 to 0.01 for R = 5 Ohms ;using 
Model transition conditions as defined by eqns(3.6b) & (3.12B) and implying to 
intermediate output of ‘ode45’ subroutine. 



FIG 3.18(a) FIG 3.18(b) 

FIG 3.18 : System Response for Change in R from 30 to 50 Ohms for D = 0. 1 ;using 
State - Space Averaged Model defined by eqns(3.4) & (3.5) . 
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factors influencing the accuracy of state space averaged model are not apparent. For 
transition where the converter remains in discontinuous conduction during the initial state, 
the transient state and the final state it is quite logical to use average model defined by 
equations 3.4 and 3.5 to study the transient, where the value of Di is taken to be (1- 

Dboundaiy)» where Dboundaiy corresponds to the value of R at final operating point. Figure 3.18 
shows the results for case 7 in Table 1 . Though the final steady state value is correct this 
method does not represent the transient correctly. 

3.4 SMALL SIGNAL APPROXIMATION FOR LINEARITY : 


If the power stage of the buck-boost converter can be linearized than bode plots and 
stability criterion can be used to determine the appropriate compensation in feedback loop 
for the desired steady state and transient response. Flow graph for linear analysis of switched 
regulators is shown in Fig. 3.19 . 

X = A0X + BOu + Ed (3.13) 

represents the converter in linearized form about a DC operating point. 

Where AO and BO are given by equations 3.2 and 3.3 or 3.4 and 3.5 and matrix E is given as 

E = (Ai-A 2 )Xo + (Bi-B 2 )V dc (3.14) 

E = (Ai - A 2 -A 3 )Xo + (Bi -B 2 -B 3 )V dc (3-15) 


I I, r, 
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SWITCHED LINEAR EQUATIONS 



FIG 3.19: Flowgraph for Small Signal Approximation for Linearity 
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depending upon the operating point lies in continous or discontinous zone of operation. 


Xo = -inv (AO) BO V dc 


V c 

3.4.1 TRANSFER FUNCTION - — AND BODE PLOT : 

d(s) 


Taking the laplace transform of equation 13 in appendix B-2 results in 


sX(s) = A0X(s) + B0u(s) + Ed(s) (3.16) 


if the input voltage V dc is fixed then u(s) = 0 

A 

Substituting this value of u(s) = 0 in equation 3.16 results in 


X(s) = (si - AO) -1 E d(s) (3.17) 


Substituting the desired values of AO and E for an operating point in continuous zone and 
equating only the first row of the matrices on LHS and RHS results in 
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(3.18) 


Vc(s) _ j_ 

d(s) “ A 


(s + ^ + 
L 


(l-D)Rr c (-RI Tn ) (l-D)R V riR 
R,L R,C R,C ^ L 


R,L + 


a -ov,v 

R, L ; 


Where A = (s+ 



q- D )Rr 1 - (1-D)’R(1-V./R,) 
R,L X R,C R,LC 


Similarly transfer function for an operating point in discontinuous operating zone can 
be derived. The transfer function shows that the power stage of the converter is of type 0. 

Fig. 3.20 to Fig. 3.23 gives the Bode plot for the operating points given in Table 2 for 
the above transfer function. 


SI. No. 

Load Resistance 

Duty Ratio 

1 . 

5 

0.1 

2. 

5 

0.5 

3. 

30 

0.1 

4. 

30 

0.5 


Table-2 

The bode plots are verified by giving a small disturbance to duty ratio d in the small signal 
model defined by equation 3.17 at the frequencies given in Table-3 Figure 3.24 shows 
the result of verification simulation at comer frequency for each case of Table 2. 
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3.4.2 TRANSFER FUNCTION ^ AND BODE PLOT : 

u(s) 


If in equation 3.16 if d(s) = 0 then 
X(s) = (si - A0)“‘ B0u(s) 


(3.19) 


substituting the values of AO and BO for an operating point in continuous zone and equating 
only the first row of RHS and LHS, results in 


Vc(s) _ (l-D)RD 
u(s) = 

Where A is given by equation 3.18. 

Similarly transfer function for an operating point in discontinuous operating zone can be 
derived. Figure 3.25 to 3.28 give the bode plot for the operating points given in Table 2 for 
the above transfer function. The bode plots are verified by giving a small disturbance to input 
voltage V dc in the small signal model defined by equation 3.19 at the frequencies given in 
Table-4 Figure 3.28 shows the result of verification simulation at comer frequency for each 


case of Table 2. 



FIG 3.20 (a): For the Operating Point R = 5 Ohms & D = 0. 1 FIG 3.21 (a): For the Operating Point R = 5 Ohms & D = 0.5 




BODE PLOTS FOR Vc (s ) / d ( s ) Transfer function 







0 15 
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BODE PLOTS FOR^c (s ) / Vdc( s ) Transfer function 
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Table 4 : Results of verification of the Small - Signal Model defined by Equation 3.19 





CHAPTER IV 


CLOSE LOOP CONTROL 


The output voltage of DC power supplies are regulated within specified tolerance 
band, also the transient and steady state response of the converter must satisfy some desired 
specifications. Section 4. 1 describes how the system pole locations determine the nature of 
transient and steady state response and how the system can be classified on the basis of pole 
location. Section 4.2 discusses how pole-placement technique can be used to design a servo 
system. Results obtained from sections 4.1 and 4.2 are used in section 4.3 to 4.5 to construct 
a servo controller for the buck-boost converter and simulate the response. 

4.1 TRANSIENT RESPONSE OF A SECOND-ORDER SYSTEM : 


Eigenvalues of matrix AO (where AO is defined by equation 3.2 or 3.4) are identical to 
the open-loop poles of the transfer function G(s). (refer to Appendix C-l). Figure 4.1 
illustrates the relationship between the location of the poles Si, S 2 , damping ratio (Z), 
damping factor (cr), undamped natural frequency (co n ) and damped frequency (o>d)- 

Si = - a + j©d 


S 2 = -o-jcod 
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When the damping factor a is positive the unit step response will settle to it’s 
constant final value. Negative damping results in a response that grows without bound with 
time and system is unstable. Zero damping results in a sustained sinusoidal oscillation and 
the system is marginally stable or unstable. 

Figure 4.2 illustrates the effect of varying the damping ratio Z from -oo to oo on the 
eigenvalue location of the second order system, vice-versa the effect of eigenvalues on the 
damping of the second-order system. 

Classification of system dynamics with respect to the value of Z is made in Table 5 


Value of Z 

System Classification 

0<Z<1 

Underdamped 

Z=1 

Critically demped 

Z>1 

Overdamped 

z=o 

Undamped 

Z<0 

Negatively damped 


Table 5 : Classification of system dynamics with respect to the value of Z 
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Step response of the system becomes oscillatory as Z is decreased, and if Z>1 the step 
response does not exhibit any overshoot and never exceeds its final value during the transient. 
Maximum overshoot of the step response of second-order system is a function of damping 
ratio Z only, while the settling time is a function of co d i.e both Z and co n . 

4.1.1 : POLE PLACEMENT 


Consider a system defined by equation 4.1 
X = AOX + BOu 


Where, 


X 

State Vector (n- vector) 

u = 

Control signal (scalar) 

A 

nxn (Constant matrix) 

B 

nxl (Constant matrix) 


(4.1) 


Regulator systems are feedback control systems that will bring non-zero states caused 
by external disturbances to the origin with sufficient speed. One approach to design regulator 
systems is to construct an asymptotically stable closed-loop system by specifying the desired 
locations for the closed loop poles. This can be accomplished by use of state feed back such 
that the control vector u = -Kx where u is unconstrained, and x is instantaneous state. 
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FIG4.3(a): State - Space Representation of Openloop System 



FIG4.3(b): State - Space Representation of Closeloop Regulator 



FIG4.4: State - Space Representation of Typel Servo System 
for a Type 0 Plant 
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Figure 4.3(a) and 4.3(b) show the open loop system and the system with state feed 
back respectively. 

The feed back gain matrix K is so decided that the system will have a desired 
characteristic equation. This design scheme is referred to as pole placement. Substituting 
u— Kx in equation 4. 1 results in 

X(t) = (A-BK)X(t) (4.2) 

and solution of equation 4.2 is given by 

X(t) = e (A ' BK)t X(0) (4.3) 

Where X(0)= initial state caused by external disturbance. If matrix K is properly 
chosen then matrix (A-BK) can be made asymptotically stable and for all X(0)^0, X(t) 
approaches 0 as t approaches infinity. 

4.2 DESIGN OF SERVO SYSTEMS : [When the plant has no Integrator] 

Servo system is a stable regulator system that will bring the error to zero given any 
initial error, between the reference and the output. The basic principle of design of type 1 
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servo system for a type 0 plant is to insert an integrator in the feed forward path between the 
error comparator and the plant as shown in Fig. 4.4. 


From Figure 4.4 


X 

= AX + Bu 

(4.4) 

Y 

= CX 

(4.5) 

u 

- (KFB)X + K,e 

(4.6) 

e 

= r-y = r-CX 

(4.7) 

Where, 

y 

Output signal (scalar) 


e = 

Output of Integrator (scalar) 


r = 

reference input signal (scalar, step function) 


C 

1 x n (Constant matrix) 



It is assumed that the transfer function of the plant has no zero at the origin that has a 
possibility of canceling the integrator being inserted. 

If r is applied at t=0 then for t>0 system dynamics can be explained by 


X(t) 


' A O' 

'X(t)' 


B 

u(t) + 

'O' 

_ e (t) 


-C 0 

|e(t)_ 


0 _ 


1 


r(t) 


(4.8) 
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(4.9) 


X(oo) 


' A O' 

'X(oo)' 


'B' 


'O' 

_ e (°°)_ 


-C 0 

_ z (°°)_ 

+ 

_0 

u(oo) + 

1 


r(oo) 


Since r is a step input 


r(co) = r(t) = r for t>0 


Therefore 


Define 


X(t)-X(«>) 


“A 0" 

rx(t)-x(co)i 


"B" 

_ e(t)-e(oo) _ 


-C 0 

1 

IT 

1 

4^ 

N «' 

<D 

1 

+ 

0 


[u(t)-u(co)] 


X(t) - X(«>) = 

X v (t) 

e(t) - e(oo) 

e v (t) 

u(t) - u(oo) 

u v (t) 


Where X v , e v , u v are extended vectors 

Then equation 4. 1 0 can be written as 


X v (t) 


'A O' 

X(t)l 

'B" 

_ e v(t)_ 


-C 0 

_e v (t)J H 

_ 0 _ 


u v (t) 


(4.10) 


(4.11) 
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Where 


u v (t)=-KX v (t)+Kie v (t) 


(4.12) 


A new (n+l) dl order vector V(t) can be defined as 


V(t) « 


X v (t) 

e v(t). 


Then equation 4.1 1 can be expressed as 


A A 

V = AV + Bu v 


(4.13) 


Where 


A 


A 0 
-C 0 



K = [-KFB-K,] 


u 


V 


A 

= KV 


The type 1 servo system results in a stable (n+l)* order regulator system that will 
bring the new error vector V(t) to zero given any initial condition V(0). 
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If the system defined by equation 4.13 is completely state controllable then by 

A 

specifying the desired characteristic equation for the system matrix K can be determined by 
pole placement (Refer Appendix C). 


4.3 OPEN LOOP DAMPING RATIO AND UNDAMPED NATURAL 
FREQUENCY 

Figures 4.5(a) and 4.5(b) show open loop system damping and undamped natural 
frequency. 

It is observed that 

(a) If the load resistance R is kept constant and duty ratio D is increased the system 
damping ratio increases and the undamped natural frequency reduces. For values of 
R>20fi (i.e. the value of R after which both continuous and discontinuous operation 
is possible) the undamped natural frequency remains constant for values of D for 
which the converter operates in discontinuous zone for given R, after which W n 
reduces as duty ratio increases. Therefore in open loop as D increases for given value 
of R system responses i.e. peak overshoot and settling time reduce inherently in open 
loop. 
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FIG 4.5 : Variation of (a) Z (b) Wn w.r.t Duty Ratio D (c) Ki (d) Ki (e) Ksuw.r.t Duty Ratio D for Z’ - 0.65 











(b) As duty ratio D is kept constant and the load resistance R is increased then for low values 
of duty ratio D the drop in system damping is not significant, but at higher values of D as 
R is increased the damping ratio decreases significantly. The undamped natural 
frequency in the continuous operation zone does not change as file load resistance R is 
increased. However if the load resistance increases to a value such that the converter 
operates in discontinuous zone then the undamped natural frequency reduces significantly 
as the load resistance R increases. 

It can be concluded that for a given duty ratio D as the load resistance increases the system 
responses (peak overshoot and settling time) increase if the converter operating point remains 
in continuous zone only, and may remain same or reduce significantly if the operating point 
transits into discontinuous operating zone. At low values of load resistance and higher values 
of duty ratio the system is inherently sufficiently damped and system responses are 
acceptable, but at low values of load resistance and low values of duty ratio the system is not 
sufficiently damped inherently. 

4.4 CLOSE-LOOP CONTROL OF THE CONVERTER 

In section 4.2 it was discussed, how state feed back can be used to construct a stable 
close-loop system for a plant represented by a linear time invariant equation 4.4, when the 

A A 

plant is of type 0 using pole-placement. The transfer function V c (s)/ d(s) in section 3.4.1 
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Dref 



FIG 4.6 : Structure of the servo controller (using pole placement) for the converter 



FIG 4.7 : Algorithm to calculate the feedback matrix KFB & Ki of the servo 
controller (using pole placement) for the converter 












shows that power stage of the converter is of type 0. The purpose of employing pole- 
placement technique to design a close-loop system is 

(a) To locate all poles of the open-loop system such that the close loop system is stable. 

(b) To locate poles of the close-loop system such that the close-loop system responses 
(i.e. peak overshoot, delay time, rise time and settling time) are within specified 
limits. 

Flow graph given in Figure shows how the state feedback pole-placement 
technique can be employed to construct a close loop servo system for the buck-boost 
converter. 

A 

To calculate K matrix , coefficients of the desired characteristic equation need to be 
known. To know these coefficients, eigenvalues of new characteristic equation 

AAA AAA 

sI-(A-BK) need to be known. Working backwards the eigenvalues of (A-BK) (i.e. new 
pole locations) are specified then the corresponding coefficients of the new characteristic 

A 

equation are determined and K matrix is calculated. Since the Open loop poles of the 
converter are already located in the left half of the s-plane, by employing pole-placement we 
can modify the system response by relocating the poles, also since an integrator has been 
inserted, thus for the new (n+1) 1 * 1 order regulator system three eigenvalues of the new 
characteristic equation need to be specified. 
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(a) 

Si 

Integrator pole location 

(b) 

S 2 

-a+j(o d 

(c) 

S 3 

-CT-jCOd 

Where 




a = 

r C0 n 


COn 



Z’ 

Desired damping ratio of closed loop system 


COn 

undamped natural frequency of open loop system 


4.5 SIMULATION OF CONVERTER WITH CLOSED-LOOP CONTROL 


4 -, £ 

Fig. ' shows the structure of the controller and it’s implementation. Fig. 3.4(a) 
shows average output voltage Vo versus duty ratio D for different values of the load 
resistance R. It is observed that in the range D=0.1 to D=0.8 for given D the output voltage 
for different values of R does not vary significantly and increases linearly with D. Therefore 
to operate the converter in the linear range of Vo-D curve such that changes in load resistance 
R do not effect the output voltage severely Vref is restricted in the range 2.5V to 35V, and 
load resistance R is restricted in the range 5Q to 50Q. 


Figs. 4.5(c) (d) and (e) show how, for the complete range of Vref and R the values of 
Ki, K 2 and Ki vary for different operating points, if desired damping ratio is chosen to be 
0.65. 
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4.5.1 CLOSE-LOOP RESPONSE FOR STEP CHANGE IN LOAD RESISTANCE R 
AND Vref 

Figure 4.6(a) and 4.6(b) illustrate the fact that for all values of load resistances and 
high values of duty ratio D the system is sufficiently damped inherently. Also it is seen that 
if Z’ the desired value of damping ratio is made 0.65 for all operating points it sufficiently 
damps the system response as seen for cases in Table 6, from Figs. 4.8 to 4.11. 


SI. No. 

Initial Condition 

Final condition 


R 


R 

Vref 

1 . 

30Q 

5V 

50fi 

5V 

2. 

10Q 

3V 

40ft 

3V 

3. 

5 Ofi 

5V 

50Q 

10V 

4. 

5Q 

15V 

5Q 

10V 


Table 6 

A 

One option is to calculate K (or KFB and Ki)corresponding to the final operating 

A 

point, other is to update the value of K according to the value of D p this can be referred to as 

A 

cycle by cycle control, another option is to put K-[l 1 1] which corresponds to the case Z >1 

A A 

or use K avera g e where K average corresponds to a feedback matrix whose elements are average 
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values of Ki, K 2 , Ki as calculated for different operating points in the range D=0 to 1 and 
R=5Q to 500. 

Fig. 4.12 shows the comparative results of using the different options 

Fig. 4.13 illustrates, how the value of KFB as calculated by the method given in Fig. 
4.6 using Z’=0.65 when used in a stand alone regulator as shown in Fig. 4.3(b) brings the 
output voltage and inductor current to reference values after they have been perturbed by an 
external disturbance. The operating point considered is R=30O and Vref=5V. 

A 

Fig.4.14 illustrates the regulator response for the above case when K=[l 1] and 

A 

K average* 

4.5.2 CLOSE-LOOP RESPONSE FOR SINUSOIDAL VARIATION IN Vref 

Fig. 4.14 shows the response for sinusoidal variation of V re f around the operating 
point defined by Vref=5V and R=30Q. 
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4.5.3 CLOSED-LOOP RESPONSE FOR SINUSOIDAL VARIATION IN V D c : 

Fig. 4. 1 5 shows the response for sinusoidal variation of input voltage Vd C around the 
operating point defined by V re f=5V and R=30Q. 


4.5.4 CLOSED-LOOP RESPONSE FOR SINUSOIDAL VARIATION IN LOAD R : 

4. !6 

Fig. 1 shows the response for sinusoidal variation in load resistance R of V re f 
around the operating point defined by Vref=5V and R=30O. 
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FIG 4.9: System Response to change in R from 10 to 40 Ohms for Vref - 3 volts 








FIG 4.11 : System Response to change in Vref from 15 to 10 volts for R = 5 Ohms 

















FIG 4.12 : System Response to change in Vref from 7.5 to 5 volts for R — 50 Ohms;using KFB — [1 1] & Ki - 
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FIG 4.13: Regulator Response for -2.5V perturbation in Capacitor voltage;using KFB = [ 1 1] & Ki 






5.5 



FIG 4.13 : Regulator Response to -2.5v perturbation in capacitor voltage at the operating point defined by 
Vref = 5 volts & R = 50 Ohms;using cycle by cycle control 
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FIG 4.14 : System Response to sinusoidal variation of Vref at 500 Hz of magnitude 2.5v at R= 30 Ohms & Vref = 20 
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FIG 4.15: System Response to Sinusoidal variation in Vdc of lv and 5000Hz around operating point Vref=5v & R=30 Ohms 
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FIG 4.16: System Response to Sinusoidal Variation in R at 50Hz of magnitude 10 Ohms ; around Vref=5v & R=30 Ohms 
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FIG 4.16: System Response to Sinusoidal Variation in R at 50Hz of magnitude lOOhms ;around Vref=5v & R=30 Ohms ;for kfb=[1 l] & ki= 





CHAPTER V 


CONCLUSION AND FUTURE WORK : 

From the modelling End the detailed analysis of the Buck-boost converter an insight 
into its characteristics has been achieved. Zones of continuous, boundary and disco ntinu ous 
operation have been defined for a given load resistance and that during the discontinuous 
mode of conduction for a given load resistance R the duration of mode-2 remains constant 
for all values of D. This information can be used to defined continuous state space averaged 
models and therefore the small signal models for any operating point easily. It is observed 
that the system satisfies the necessary and sufficient conditions of arbitrary pole placement 
and state feed back controller has been designed. 

However, the effort to prove the fact that the state space averaged models can 
represent the actual converter for all kinds of perturbation has succeeded partially. For the 
transition of the kind DCM — DCM — DCM, DCM — CM — DCM, CM — DCM — DCM during 
the initial state, the transient state and the final state additional factors governing the model 
transition conditions need to be examined. One possibility is to examine the nature of 
transients in the diffrential equation model with particular attention being given to variation 
of Di. If this variation can be modelled then probably the information gained can be used in 
the state space averaged models to represent all kinds of perturbations. 



Deciding the desired damping ratio to achieve optimum system response for change 
in V re f, Vac, R is difficult. Fuzzy decision making can be applied to decide close loop 
damping ratio by fuzzyfying the open loop response. It is observed that values of Ki, K*,, Ki 
as desired by the algorithm given in Fig. f. gives satisfactory results. However optimization 
of the response can be done by using fuzzy gain sched uling 
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State Vector . A vector that determines uniquely the system state x(t) for any time t> 
t 0 once the state at t=t 0 is given and the input u(t) for t> t 0 is specified. 

State Space : The n-dimensional space whose co-ordinate axes consist of the xi , X2 
. . x n axis. Any state can be represented by a point in the state -space. 

State Space Equations : Assume that a multiple input -output system involves n 
integers, also there are r inputs Ui (t) , U2 (t), ...,u r (t) and m outputs yi (t) , y 2 (t), 
. . . ,y m (t) 

Let the n outputs of the integers be defined by state variables xi (t) , X2 (t), ...,x n (t). The 
system may then be described by 

Xi (t) = f, (Xi , X r ; Ui , u 2 , . . .,u r ;t) 

X2 (t) = f 2 (Xi , X 2 >,Xr j Uj , U2 , ..^Ur 9 t) 

X n (t) = f n (X! ,X 2 ,...,X r ;U] ,U 2 , ...,U r ;t) 
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and the system outputs given by 


y,(t) = g| (xi , x 2 ,...,x r ; ui , u 2 , ...,u r ;t) 


y 2 (t) = g,(xi ,x 2 ,...,x r ;ui,u 2 


..,u r 


;t) 


y m (t) = g m (xi , x 2 ,...,x r ; U) , u 2 , ...,u r ;t) 


Xl(t) 


X(t) = 


x 2 (t) 


Xn(t) 



'fi(x„x 2) . 

...X^Uj,^,.. 

..,U r ,t) 

f(x,u,t) = 

f 2 (x,,x 2 , 

..,x„;uj,u 2 ,.. 

...,u r ;t) 


_ f n( X P X 2>- 

....,x n ;u p u 2 u r ;t)J 


Y(t) = 


y,(t) 

y 2 (0 


y n (t) 


u(t) = 


U,(t) 

u 2 (t) 


Un(t) 


Then, X(t) = f(x,u,t) 


g(x, u, t) = 


g,(x 1 ,x 2 ,...,x n> u 1 ,u 2 , ,u r ,t) 

g 2 (x 1 ,x 2 ,...,x n ;u„u 2 ,u r ;t) 

g n (x 15 X 2 , ,X n ,U],U 2 ,....,U r ,t) 


( 1 ) 


Y(t) = g(x, u,t) 


( 2 ) 
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Eq. 1 is the input Equation, and Eq. 2 is the output Equation. If vector function f and/ or g 
involve time t explicitily then system is called a time -varying system and if Eq. 1 and 2 are 
linearised about the operating point then we have the following linearized eqn. 

X(t) = A(t)X(t) + B(t)u(t) 

Y(t) = C(t)X(t) + D(t)u(t) 

Where A(t) is called the state matrix 
B(t) is called the input matrix 
C(t) is called output matrix 
D(t) is called direct transmission matrix 
If functions f and g do not involve t explicitly then system is called time variant system and is 
represented by 

X(t) = AX(t) + Bu(t) (5) 

Y(t) = CX(t) + Du(t) (6) 

Fig. A-l gives the block diagram representation of equations. 


(3) 

(4) 



APPENDIX B 


B-l CONCEPTS OF STATE-SPACE AVERAGING : 

Let us consider the two method conduction of the converter: The system can now be 
represented by 

X(t) = A,X(t) + B,u(t) for n T s < t< (n+D)T s 
X(t) = A,X(t) + B 2 u(t) for (n+D) T s < t< (n+l)T s 

for any cycle n. 

Where T s =l/f s 

Eq. 1 cannot continuously follow Eq. 2 . Since X [ (n+D)T s -s] as described by Eq. 1 cannot 

equal X[ (n+D)T s ] as described by Eq. 2, no matter how closely a approaches zero. But if 
the switching components of the X(t) waveform are excluded continous description of a 
discontinuous system can be approached. Now the time-interval of Eq. 1 can include 
(n+D)T s and that of Eq. 2 can include (n+l)T. . If <he switching frequency is high enough 


( 1 ) 

( 2 ) 
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relative to the circuit natural frequency and the signal frequency then switching period is 
short enough to approximate X at t = n T s and at t = (n+D)T s as follows 


X(nT„) = 


X[(n + D)TJ - X[nT„1 
DT S 


X(n + D)T S 


X[(n+1)TJ-X[(n+D)TJ 

0-D)T s 


( 3 ) 

( 4 ) 


Since now closed time intervals are permitted, we combine Eq. (1 and 3) & (2 and 4), to 
eliminate X(nT s ) , X[(n + D)]T S 


X(n + D)T S = X(nT s ) + DT S [ A,X(nT s ) + Bi u(nT s )] (5) 

X(n + 1)T S = X[(n+D)T S )] +(1+ D)T S [ A 2 X(n+D)T s + B 2 u(n+D)T s ] (6) 

Combining Eq. 5&6 to eliminate X [ (n+D)T s ] we get 

X[(n+1)T S ] = X (nT s ) + DT S [ AiX(nT s ) + B lU (nT s ) + 

(1-D) T s A 2 [ X(nT s ) +DT S [A,X((nT s ) +Bju (nT s )]]+ (1-D) T s B 2 u (n+D)T s (7) 

„ T u[(n + D)TJ - u(nT s ) 

Now u(nT.) = — 

DT S 



Therefore 


X[(n + 1)T S ] - X(nT s ) 

— ~ - [A|D+ (l-D)A 2 ]X(nT s ) + [B,D+ (1-D)B 2 ]X(nT s ) 

+ [(l-D)D T s [ A 2 [A]X(nT s )+B, u(nT s ) ]]+ [B 2 u (nT s ) ] (8) 

If T s is very small LHS of Eq. 8 approximates X(nT s ) and T s terms of RHS being negligible 
results in 

X (nT s ) = [AjD+ (l-D)A 2 ]X(nT s ) + [B,D+ (1-D)B 2 ]u(nT s ) (9) 

Therefore 

A0 = A,D + A 2 (1-D) 

B0 = B,D + B 2 (1-D) 

B-2 LINEARIZATION OF STATE -SPACE AVERAGED MODELS : 

In the linear state space averaged Eqn,. if the D.C. terms of X,u, D are separated from 
signal frequency A-C terms. It is assumed that ac amplitudes are small enough so that the 
product of any two ac terms is negligible. 
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Substituting X = X 0 + X 

A 

U = Uo + u 

A 

D = D 0 + d 

AAA 

(where X, u, d identify the signal frequency ac terms) in Eq. 9 results in DC Eqn. 

0 = [A,D+ (1-D)A 2 ]Xo + [B&+ (I-DJB 2 K (10) 

and ac Eq. 

X= [A, EH- (1-D)A 2 ] X + [B,D+ B 2 (l-D) ] u + [(A,-A 2 )]X 0 + [(Bi-B 2 ) uj d (11) 

Defining E = (Ai-A 2 )X 0 + [(Bi-B 2 )]u 0 

Eq. 10 and 1 1 can be written as 


0 = AO X 0 + BOUo 

(12) 

A A A 

X = AOX + BO u +Ed 

(13) 
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Adding Eq. 12 & 13 results in a linearized Eqn. 

X = AOX + BO u + Ed 

Valid for an operating point in the continous operating zone. But in discontinous operating 
zone AO and BO are defined by Eq. 3.4 and 3.5 and matrix E is defined as 

E = (A r A 2 - A 3 ) Xo + (Bi - B 2 -B 3 )Uo (14) 


Where X 0 = -inv (AO) BO Vd C 


B-3 SMALL SIGNAL MODEL OF PULSE WIDTH MODULATOR : 


In the direct duty ratio pulse-width modulator, the control voltage v c (t) , which is 
output of the error amplifier is compared with a repetitive waveform v r (t), which establishes 
the switching frequency f s , as in Fig. B-J~ The control voltage v c (t) consists of a dc 
component and a small ac perturbation component 


v c (t) = v c + v c 


(15) 
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Vc(t) 



Graphical Representation of functioning of Pulse Width Modulator 


where v c (t) is in a range between 0 andv r , as shown in Fig. B-l. v c (t) is a sinusoidal ac 


perturbation in the control voltage at a frequency to, where a is much smaller than the 

switching frequency co s (— 2IIf s ). The ac perturbation in the control voltage can be expressed 

as 

v c (t) = a sin (<nt-<J>) 

(16) 

by means of an amplitude a and an arbitrary phase angle <j). 


In fig. B-2 , the instantaneous switching duty ratio d(t) is as follows : 


d(t) =1 .0 if v c (t) > v r (t) 

(17) 

= 0 if v c (t) < v r (t) 

(18) 


d(t) can be expressed in terms of the Fourier series as 

d(t) =— + -£- sin(cot - <(>) + other high frequency component ( 1 9) 

v r Vr 

The higher frequency components in the output voltage v 0 due to the high frequency 
components in d(t) are eliminated because of the low pass filter at the output of the converter, 
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therefore the high frequency components in Eq. 13 can be ignored. In terms of its dc value 
and ac perturbation : 


d(t) = D+ d(t) 


( 20 ) 


Comparing Eq. 19 & 20 



( 21 ) 


d(t) = — sin(oot — (|) ) 
V, 


( 22 ) 


From Eq. 16 & 22, the small signal transfer function T m (s) of the Pulse Width Modulator is 


T m (s) = 


d(s) 


Vc(s) 


1 
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appendix c 


C-l : CORRELATION BETWEEN TRANSFER FUNCTION AND 
STATE-SPACE EQUATIONS : 


Equations 5 and 6 in Appendix A represents a system in state-space. Since transfer function 
is defined as the ratio of the Laplace transform of the output to the Laplace transform of the 
input when the initial conditions are zero, then taking X(0) as zero in Laplace transform of 
the equation (5) in appendix A results in: 

sX(s) - X(0) = AX(s) + Bu(s) 
sX(s) - AX(s) = Bu(s) 

(sI-A) (Xs) = Bu(s) 

Multiplying both sides by (sI-A)' 1 

X(s) = (sI-A)' 1 Bu(s) 

Substituting this in equation 6 Appendix A results in 
Y(s) = [C(sl-A)' 1 B+D] u(s) 

And the transfer function G(s)=C(sI-A)"' B+D (1) 
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Right hand side of equation 1 involves (si- A)'. Hence G(s) can be written as 
G(s) = Q(s)/ jsI-Al 

Where Q(s) is a polynomial in s. Therefore |sl - A| is equal to the characteristic polynomial 
of G(s) and the eigenvalues of A are identical to poles of G(s). 

C-2 : SOLUTIONS OF LINEAR TIME INVARIANT VECTOR MATRIX 
DIFFERENTIAL EQUATIONS : 


Consider the following vector matrix differential equation. 

X(t) = AX(t) + Bu(t) 

Let X(o) = X 0 

Where X is an n-vector 

u is a r-vector 
A is nxn constant matrix 
B is nxr constant matrix 
Equation 2 can be written as 

X (t)-AX(t) = Bu(t) 
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Multiplying both sides of this equation by e At results in 


e A, [X(t)~AX{t)J = i-[e' A 'X(t)3 = e - A, Bu(t) 


integrating the above equation between 0 and t 


e~"X(t) = X„ + Je"Bu(t)d(t) 

0 


X(t) — e A *X 0 + Je Ait ?) Bu(i)di 
0 


( 3 ) 


If initial time is to instead of o then 


X(t) = e IA(I '•*X # + Je A ‘ , - t, Bu(T)(h 
0 


Where X* = X(to) 

The response X(t) is a sum of the motions due to the initial condition and those due to 
the forcing function. Motion due to forcing function depend on B. It is possible to choose B 
such that any particular motion e x,t cannot be excited for any input. 
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C-3: COMPLETE STATE CONTROLLABILITY OF CONTINUOUS 
TIME SYSTEMS 

Solutions of linear time-invariant vector matrix differential equation defined by 
equation 2 is given by equation 3. Note that using the^^ster’ s m terpolatio^toula for a 
minimal polynomial involving only distinct roots e' At can be written as ^ ' 

e ' AT *Z a k( T ) AK (4) 

k 0 

System described ■ equation 2 is said to be state controllable at t=tO if it is possible 
to construct an unconstrained contract signal that will transfer an initial state to any final state 
in a finite time interval tO<t <t]. If every state is controllable them the system is said to be 
completely state controllable. 

CONDITION OF COMPLETE STATE CONTROLLABILITY : 

It is assumed that the final state is the origin of the state space & that the initial time 
t0=0. Applying the definition of complete state controllability in equation 3 results in : 

«i 

X(t,) = e At, X 0 + Je A(tl_t) Bu(T)dt = 0 (5) 
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(6) 


‘i 

X(0) = -jV At Bu(x)dx 

0 

Substituting the value of e' At as given by equation (4) 

n-1 h 

X(0) = -^A k B fa K (x)u(x)dx 

k=0 o 

t 

Put ja K (x)u(x)dx = P K 

0 

X(0) = -Ea k B|3 k 

lc=0 

Po 
p, 


If the system is completely state controllable then given any initial state X(0), equation 6 
must be satisfied. This requires that the (rank of th e nxn matri ^ 1 consisting of(iiJjn eari}( 
independent colmnn vectors) [B AB A n+1 B] (called as the controllability matrix) be n. 

C-4 NECESSARY & SUFFICIENT CONDITIONS FOR ARBITARY 
POLE PLACEMENT.: 

Consider the control system defined by equation.2. We assume that the magnitude of 
the control signal u is unbounded. If the control signal u is chosen as : u=-Kx, where K is the 



X(0) = -[B AB A n_1 B] 
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state feedback gain matrix (lxn matrix), then the system becomes a closed loop control 
system as shown in Fig. 4.3(a) and the solution to equation.2 becomes 

X(t)=e (A ' BK)t X(0). 

Note that the eignvalues of matrix A-BK (which we denote pi , p 2 , u n ) are the desired 

closed loop poles. If the system is not completely state controllable, then there are 
eigenvalues of matrix A-BK that cannot be controlled by state feedback. Suppose the system 
represented by equation 2 is not completely state controllable. Then the rank of the 
controllable matrix is less than n, 

or rank [B|AB| A n_1 B] = q < n 


This means that there are q linearly independant column vectors in the controllibility matrix, 

such q linearly independent column vectors are defined as fi, £ 2 , f q . Also, let n-q 

additional vectors V q +i, V q+ 2 V n be such that 


Then 


P = [f,f 2 f q V q+1 V q+2 VJ is of rank n. 


A = P~'AP = 




B = P _1 B = 



Define K = KP = [k, K 2 ] 
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Then 


si - A + BK| = |P"' (si - A + BK)P| 

= | si - P _1 AP + P _1 BKP| = si - A+ B K 

A A A 

sI-A + BK|= sI-A+BK 



slq An+BuK] A 12 + BjjKj 


|sl q A u + BjjKjj . |sl n _ q A 22 | — 0 

Where I q is a q-dimensional identity matrix hq is an (n-q) dimensional identity matrix. 

Notice that the eigenvalues of A 22 do not depend on K. Thus if the system is not completely 
state controllable, then there are eigenvalues of matrix A that cannot be arbitrarily placed, 
therefore, to place the eigenvalues of matrix A-BK arbitrarily, the system must be completely 
state controllable (necessary condition) 

Next, if the system is completely state controllable [means that matrix M given by 
equation 9 below has rank n or has inverse], then all eigenvalues of matrix A can be arbitanly 
placed (sufficient conditions). In proving the sufficient condition, it is convenient to 
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transform the state equation given by equation (2) into controllable Canonical form. Define a 
transformation matrix T by 


T=MW 


( 8 ) 


Where M is the controllability matrix 


M = [B|AB .|A n -'B] 


( 9 ) 


and 


‘n-l a n-2 

a, 1 

n-2 a n-3 

a, 0 


0 

o o 

1 

1 

00 


Where the &\ ,s are co-efficients of the characteristic polynomial. 


|sl- At = S" + ai S n + a, S"-‘+ .+a n _,S + a n 


Define a new state vector X by 


A 


X =TX 



If the rank of the controllability mark M is n (meaning that the system is controllable), then 
the inverse of matrix T exists and equation (8) can be modified to 


X = T“ 1 ATX+T~ 1 Bu 


r'AT 


0 1 0 

0 0 1 

0 0 0 

— ^n-1 _ ^n-2 


0 

0 

1 

a, 


0 

0 


o 


1 


( 10 ) 


Equation (10) represents the system in controllable canonical form. Thus, given a 
state equation. Equation (2), it can be transformed into the controllable canonical form if the 
system is completely state controllable and if we transform the state vector X into state vector 

A 

X by use of transformation matrix T. 
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Let the set of the desired eigenvalues be m, u 2 , 


,u n . Then the desired 


characteristic equation becomes 


(S-ui)(S-u 2 ) (S-u n ) - SVaiS”' 1 + +a n .iS+an (11) 

Let, 

K = KT = [5 n 5 n -, 8,) (12) 

A A A 

When u = - KX = -KTX is used to control the system given by equation (2) the system 
equation becomes 

X = T" 1 ATX-r 1 BKTX 

The characteristic equation is 


sI-r 1 AT+r’BKT| = 0 


This characteristic equation is the same as the characteristic equation for the system, defined 
by equation (2), when U=-KX is used as the control signal. This can be seen as follows. 

Since 

X = AX + Bu = (A - BK)X 


117 



The characteristic equation for this system is 


|sI-A + BK| = |r' (si - A + BK)T| = | si - T~‘AT + T‘ l BKT| 


Simplifying the characteristic equation of the system in the controllable canonical form 


sI-r'AT + r'BKTl 



0 

1 

0 ' 


' 0 " 


si — 

0 

0 

1 


0 

[5A-! 5,] 


" a n 

a n -l — 

“ a l_ 


_ 1 _ 



s -1 0 

0 s 0 

0 

a „ +5 n a r,-l + 5 n-l S + ^+SJ 


s n + (a t +5,)s n - 1 + .+(a n _, +5 n _,)s+(a n + 8 n ) = 0 


This is the characteristic equation for the system with state feedback. So it must be 
equal to equation (1 1) the desired characteristic equation. By equating the coefficients of like 

powers of s. 
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a, +8j = <*! 

®2 ^2 = ®2 

i 

i 

a„ +5 n = a n 

n n n 

Solving the preceding equations for the 8i s and substituting them into equation (12), results 
in 

K = KT _1 = [8 n 8 n _, 8JT" 1 

Thus if the system is completely state controllable, all eigenvalues can be arbitarily placed by 
choosing matrix K according to equation. We have thus proved that die necessary and 
sufficient condition for arbitary pole placement is that the system be completely state 
controllable. 

C-5 : 

"A 0] 

If matrix P is defined as and has rank (n+1) then the system defined by 

v U 

equation 4.13 is completely state controllable as shown : 

IfM is defined as 

M = [BAB A"’ 1 B] 


119 



and if system is defined by 


X = AX + Bu 


is completely state controllable and rank of matrix M is n then rank of 


M O' 
0 1 


is n+1 


Consider 


1 

> 

w, 

"M 

o' 


' AM 

B' 

L-C oj 

0 

1 


- CM 

0 _ 


(13) 


Since matrix P is of rank n+1 the left side of equation (13) is of rank (n+1) and 


RHS of equation (13) is also of rank (n+1), since, 


' AM 

B" 


-CM 

L 

0 



A[B AB~ A n_1 B] 

-C[B AB A n_1 B] 


B 

0 
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and if system is defined by 


X = AX + Bu 


is completely state controllable and rank of matrix M is n then rank of 


M 0 


Consider 


' A B" 

'M 

0 " 


' AM 

B“ 

-C 0 

0 

1 


-CM 

0 _ 


Since matrix P is of rank n+1 the left side of equation (13) is of rank (n+1) and 


RHS of equation (13) is also of rank (n+1), since, 


' AM B' 


'A[BAB" A n_1 B] B" 

-CM 0 


- C[B AB A n_1 B] 0_ 


(13) 
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Where rank of 


' AM B 


ABA 2 B" A n B B 

-CM 0_ 


- CB-CAB CA n_, B 0 


[ABA 2 B A n B.B] 


A A A A A A A 

[BABA 2 B A n B] is (n+-l) 


Thus the system defined by equation 4.1 3 is completely state controllable and solution 
to it can be obtained by pole placement approach. 


The state error equation can be obtained by- substituting u v = -KV in equation 4.13. 


A AAA 

V = (A-BK)V 


If the desired eigenvalues of matrix (A- BK) or the desired closed-loop poles are 

specified as ui, U 2 , u n +t then feedback gain matrix KFB and the integral gain constant Ki 

can be determined. 
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a i> a 2 a p are the co-efficients of the characteristic equation si- A and on, 

a 2 j ,otp are co-efficients of the desired characteristic equation formed by defining the 

desired eigenvalues ui, U 2 ,u p . 

Then, 

AAA A A 

M = [B AB A n_1 B] 

-- a, r 

-a, 1 0 

0 0 0 

0 0 0 

T = MW 

K = [K, ,K 2 , K p _, - K,] = [a p - a p a p _, - a p ., a, - a JT 1 

=[KFB - Ki] 

Where KFB = [KiK 2 K M ] 
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